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Simulations of the kinetic friction due to a layer of adsorbed molecules between two crystalline surfaces are presented. 
The adsorbed layer naturally produces friction that is consistent with Amontons' laws and insensitive to parameters 
that are not controlled in experiments. The kinetic friction rises logarithmically with velocity as in many experimental 
systems. Variations with potential parameters and temperature follow variations in the static friction. This correlation 
is understood through analogy with the Tomlinson model and the trends are explained with a hard-sphere picture. 
Keywords: kinetic friction, Amontons' laws, surface layers, molecular dynamics, simulations 
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1. Introduction 



One of the major challenges for researchers in the 
field of nanotribology is to uncover the molecular origins 
of macroscopic behavior. For example, why do many 
systems obey Amontons' laws that friction is propor- 
tional to normal load L and independent of the appar- 
ent geometric area of the surfaces Aapp? Unfortunately 
many simple molecular scale models give results that 
are qualitatively inconsistent with macroscopic mea- 
surements [||,i||||,|,|0,|ll| . Indeed they suggest 
that static friction should almost never be observed 
[|§M|MMl|! A possible origin for such discrepan- 
cies [^12^]lJ] is that the models consider contacts 
between two bare surfaces. Real surfaces almost al- 
ways have an intervening layer of material called " third 
bodies" p5| , p^ , p7[ . This layer may be dust, grit, wear 
debris, or the few angstroms of hydrocarbons and water 
that adsorb rapidly on surfaces that are exposed to air. 

In a recent paper He et al. explored the effect of 
adsorbed layers on the static friction Fg. They showed 
that such layers naturally lead to a non-vanishing Fg 
that exhibits much of the behavior seen in macroscopic 
experiments. For example the calculated value of Fs is 
insensitive to many parameters that are not controlled 
in experiments, including crystallographic alignment, 
the thickness of the adsorbed layer, and the precise size 
of the molecules. More importantly, the simulation re- 
sults provided a molecular basis for Bowden and Tabor's 



local pressure P: 



T = To + aP 



(1.1) 



phenomenological explanation of Amontons' laws |18 

Bowden and Tabor had noted that Aapp is usually 
much larger than the area of intimate molecular contact 
^reai- They pointed out that both Amontons' laws and 
many exceptions to them could be explained ?/the local 
shear stress r in the contacts increased linearly with the 



Summing over A^cai and dividing by load gives a friction 
coefficient 



H = F/L = a + To/P. 



(1.2) 



This is independent of load if P is constant or tq/ P <C 
a. Later work showed that both elastic||l^,0 and 
plastic [|l^ models of surfaces with many contacts yield 
a constant P at high loads. At low loads, P decreases, 
explaining why many systems show an increase in /i in 
this hmit |l|,|ll- 

He et al.'s work |l^ showed that adsorbed layers lead 
to a static friction that obeys Eq. However, most 
experiments measure the kinetic friction. Although ex- 
perimental values for the two quantities are often cor- 
related, it is not obvious that they should be. Static 
friction is the threshold force that must be overcome 
before one solid slides over another. It arises because 
the two bodies have locked together into a local energy 
minimum, and Fg is the force needed to break them 
free from this minimum. Kinetic friction Fk(v) is the 
force needed to maintain sliding at a given velocity v. 
As such it is proportional to the work that must be 
done due to energy dissipation during sliding. Thus 
kinetic and static friction involve different phenomena 
and they have very different values in many simple mod- 
els @,|,il,|,|0Jii,0|l). 

In this paper we consider the effect of adsorbed lay- 
ers on kinetic fr ictio n. We first show that the linear 
relation in Eq. 1.1 is obeyed for each velocity. The 



intercept tq is relatively independent of velocity below 
about 1 m/s. The slope a shows a logarithmic depen- 
dence on velocity over more than two decades in ve- 
locity. This same functional dependence is observed in 
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experiments, and used in rate-state models of friction 
[ p2|p^ , |2^ . We next examine the effect of interaction 
potentials, temperature, and crystallographic alignment 
on kinetic friction. The variations with these parame- 
ters are consistent with a simple hard sphere model for 
lateral forces. We conclude by considering the strong 
correlation between calculated values of Fk and Fg . A 
qualitative explanation for this correlation is suggested 
based on results for the Tomlinson model and this 
analogy is confirmed by detailed analysis of the dynam- 
ics of individual molecules. 

The following section describes the model we use for 
the surfaces and the intervening molecular layer. The 
next section presents results, and the final section pro- 
vides a summary and conclusions. 



2. Computational Methodology 

Since our goal is to establish general mechanisms, we 
use relatively simple interaction potentials in our sim- 
ulations. These allow us to quickly span a wide range 
of sliding velocities, system sizes, geometries and inter- 
actions. They also allow treatment of longer time and 
length scales than more detailed potentials. 

We find similar behavior for both crystalline and 
disordered walls For the results described 

below, each solid contains discrete atoms forming the 
(111) surface of an fee crystal. The coordinate system 
is chosen so that the surfaces lie in the x — y plane, and 
periodic boundary conditions are applied in these direc- 
tions to minimize edge effects. Note that the periodic 
boundary conditions also prevent us from considering 
perfectly incommensurate systems. However, we find 
that surfaces exhibit incommensurate behavior as soon 
as the common period becomes longer than a few lattice 
constants 

The two opposing surfaces have approximately the 
same nearest- neighbor spacing (i„„, and the effect of 
commensurability is explored by rotating the top sur- 
face by an angle 9 about the z-axis. A discrete set of 
angles is studied that allow both surfaces to retain their 
hexagonal symmetry. Except where noted, the bottom 
wall has dnn — 1.2ct. The value of d„„ for the top wall 
is adjusted slightly (< 2.1%) in order to satisfy the pe- 
riodic boundary conditions. No correlation between the 
size of this adjustment and the calculated friction was 
seen at the level of the statistical fluctuations. 

Atoms in the bottom solid are fixed in place, and 
atoms in the top solid translate together as a rigid body 
of mass M. Previous work shows that including elas- 
tic deformations of the solids produces a slightly lower 
static friction Jl2| , p"3|j25|] , and we find a similar reduction 
in the kinetic friction. A constant downward pressure is 
applied to the top solid, and its height varies in response 
to this external pressure and molecular interactions. A 
constant lateral force or velocity is applied along the 



nominal sliding direction, usually x. The wall is free 
to move in the other lateral direction, i.e. the external 
force is zero in this direction. The kinetic friction is de- 
termined by imposing a constant lateral velocity u, and 
calculating the average opposing force. This average is 
performed over a displacement of at least 25d„„ in order 
to reduce statistical fluctuations. The static friction is 
obtained by applying a constant lateral force, and de- 
termining the value needed to initiate sliding fl^ . 

The layer of molecules between the solids is described 
with a bead-spring model. Previous studies have shown 
that this model yields realistic dynamics for polymer 
melts [^6|, and shown how to map between it and de- 
tailed chemical models of polymers [p7|j2^ . These de- 
tailed models take at least an order of magnitude more 
computer time. 

Each molecule contains n spherical monomers of 
mass m. We considered simple spherical molecules 
n = 1, and short chains that model airborn hydrocar- 
bons, n = 3 and 6. All monomers interact through a 
truncated Lennard- Jones potential pof : 

VLj{r) = 4e [((7/r)i2 - {a/rf] for r < r" , (2.1) 

where r is the distance between molecules and the po- 
tential is zero for r > r'^. Adjacent monomers along a 
chain are coupled by an additional potential that pre- 
vents chains from crossing or breaking: 



VcH{r) = -{l/2)kRl\n[l - {r/Rof 



(2.2) 



where Ro = l.Scr and k — SOe/cr^ ||2^. The parameters 
e, CT, and t^j = (a^/me)^/^ are characteristic energy, 
length, and time scales, respectively. These character- 
istic scales are used to normalize other quantities. Val- 
ues that would be representative of hydrocarbons are 
e ~ 30meV, a ^ 0.5nm and t^j ^ Bps The miit of 
pressure, ecr~^, corresponds to about 40MPa. 

Wall atoms interact with monomers through a 



Lennard-Jones potential with a different cutoff, 
and different energy and length scales, twf and CTu,/, 
respectively. This allows us to increase or decrease the 
amount of adhesion, and the effective surface corruga- 
tion of the walls. In most cases we considered purely 
repulsive interactions: r'^ — 2^1^ a and r^^ — i^/^a^f. 
In this limit, atoms on opposing walls are far enough 
apart that they do not interact. As shown below, the 
main effect of increasing the interaction range is to add 
an extra adhesive pressure between the walls. 

All of our runs are done in a constant temperature 
ensemble. We have performed extensive tests on the 
effects of thermostats 30 , 3l| , |32| . In the appropriate 
limits, thermostats prevent a gradual rise in tempera- 
ture due to frictional dissipation without having a direct 
influence on the friction force. To ensure that one is in 
these limits, the thermostat must be applied normal to 
the sliding direction and be sufficiently weak, and the 
velocity must be low (< Q.la/tLj). At high velocities, 
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Figure 1. Kinetic shear stress ti^ as a function of pressure for the 
indicated velocities. Velocities are in units of cr/ti^j ~ lOOm/s, 
and the unit of stress and pressure is ea~^ ~ 40MPa. For these 



runs 9 



30° 



6, dn 



1.2a, and 



the number of monomers is equal to the number of atoms on 
the surface of the bottom wall. The data for each velocity lies 
on a straight line. We confirmed that linearity extended up to 
P/e(T~^ = 100 for the highest velocities. 

the structure and other properties of the film may be af- 
fected [ ^3p^ . In the following simulations, a Langevin 
thermostat with damping rate 0.4t£j is coupled to the 
equations of motion in both directions normal to the 
sliding direction 
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The equations of motion are integrated using a fifth- 
order Gear predictor-corrector algorithm with a time 
step dt = 0.005a/tLj- Unless otherwise noted, the area 
of a periodic unit cell is ^ = 720cr^ and the mass of the 
wall is M/m = 0.8A/a^. Previous studies show that 



finite-size effects are very small at this system size [25 



Most simulations are done at temperature T — O.Tfcs/e, 
where ks is Boltzmann's constant. 



are plotted against the logarithm of velocity in Figure 
|. The value of a increases as the logarithm of veloc- 
ity over the entire range of Fig. ^ ^ a factor of 200 in 
V. In contrast, the value of tq becomes independent of 
velocity for v < O.OIct/^lj ~ Im/s. It is negative be- 
cause these simulations used a purely repulsive poten- 
tial r^/cr = 2^/^. We show later that including adhesive 
interactions leads to positive values of tq. 

A logarithmic variation of friction with velocity is 
also seen in many experimental systems [^,^ 24 . 
Measured friction coefficients fit a "rate-state" equa- 
tion: 

/i = a + To/P = MO + Aln{v/v°) + B ln(e/e"), (3.1) 

where is a variable that describes the state of the 
system, and fiQ is the friction at the reference velocity 
v'-' and state 0". In most cases Q is treated as a phe- 
nomcnological parameter whose evolution is described 
by: 



de/dt = 1 - &v/Dc 



(3.2) 



where Dc has units of length. Dieterich and Kilgore 
have shown that changes in O are proportional changes 
in the true area of contact between rough surfaces and 
that Dc corresponds to the typical diameter of contacts 

. . . 

Our simulations give the dependence of friction on 
V at fixed contact area ||3^ and thus the slope of Fig. 
|(a) gives A = 0.00110 ± 0.00005. This is smaller than 
typical experimental values for glass or rocks where A 
is .005 to .015. However, one might expect A to in- 
crease with a. While our values of a are comparable 
to values for fiat crystalline surfaces like mica|3^,^ or 
MoS2[|l^, they are much smaller than values for rocks. 
The ratio A/a is about .05 for our system, and varies 
from .008 to .025 for rocks. We do not know of any 
measurements of A for mica or M0S2 , but it would be 
interesting to see how results for these systems compare 
to our simulations. 



3. Results 

3.1. Variation with pressure and velocity 

Figure |l| shows the pressure dependence of the kinetic 
friction at the indicated velocities for an incommensu- 
rate system, 9 = 30°. At each velocity the shear stress 
T/c rises linearly with pressure over the entire range of 
pressures. Using realistic values of e and a, the highest 
pressures correspond to about 1.5GPa. The linear rela- 
tion continues up to P = lOOecr"'^ in the cases we have 
tested. However, the larger forces at these pressures re- 
quire a reduction in the integration time step, making 
the calculations much slower. 

Values of a and tq were obtained from linear fits to 
data like that shown in Fig. 0. In each case the cor- 
relation coefficient was greater than .999. The results 



3.2. Variation of fi with interaction potential, T and 
geometry 

A variety of interaction potentials and geometries 
have been considered to determine which factors infiu- 
ence a and tq and which leave them unchanged. As 
shown in Fig. |^, doubling the interaction between wall 
atoms and fiuid monomers has almost no effect on the 
friction. In contrast, the friction increases rapidly with 
the ratio of d„„ to a^f. For dnn = l-2cr we find a rises 
from 0.012 to 0.034 as cr^/ drops from 1.2 to 0.9. The 
value of To increases in magnitude from -0.07 to -0.16. 
Although not shown, increasing d„„ at fixed a^f also 
increases a. 

These trends in a can be understood from a simple 
geometrical picture. At the pressures of interest, the 
interactions are dominated by the repulsive r~^^ term 
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Figure 2. Plot of (a) a and (b) tq against the logarithm of the 
velocity for the system of Fig. |l| 
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Figure 3. Plot of shear stress against pressure at d = 0.005a/tLj 
for the indicated values of a^f and e^j/- Other parameters are 
the same as for Fig. 0. 



in the Lennard- Jones potential. Monomers and wall 
atoms can be thought of as hard spheres with an effec- 
tive diameter D that is determined by the separation 
r at which the Lennard- Jones force balances the force 
from the applied pressure. This picture gives a diame- 
ter D cx awfi^wf /o'wfP)~^^^^ , that varies rapidly with 
(Twf but very little with e^f or P. As illustrated in 
Fig. |[ monomers can penetrate deep between surface 
atoms when D/dnn is small. This increases the friction 
because monomers must be lifted up the ramp formed 
by the surface of closest approach in order to translate. 
The required lateral force is the normal load times the 
slope of the surface which increases as D/dnn decreases. 
Thus the hard-sphere model explains the linear depen- 
dence of T on pressure as well as the insensitivity to Cuif 



I nn I 



Figure 4. Sketch of surface of closest approach for D/d„„ = 0.7 
(solid line) and 1.3 (dashed line). The lateral force required to 
lift monomers up the ramps that these surfaces represent is given 
by the slope times the normal force. 



and the variation with a^f. 

The simulations just described were for strictly repul- 



sive potentials 



'wf 



i^^^o'wf), and one may wonder 



whether the hard-sphere picture is more generally ap- 



plicable. We have done simulations with r'. 



wf 



2.2(7. 



wf 



to test this. When direct wall/wall interactions are 
neglected, we find that the kinetic friction is fit by 
a = 0.0211 ±0.0003 and tq = -h0.007± .008 over a pres- 
sure range from to 36ecr~'^ for the parameters of Fig. 
|l|. This value of a is only 8% smaller than for a purely 
repulsive potential. In contrast, tq increases substan- 
tially because the attraction between monomers and 
wall atoms increases the effective pressure holding them 
together. Adding direct wall/wall interactions leads to 
a slightly larger adhesive pressure, and a slightly larger 
value of To . 

The effect of temperature on a is also fairly weak. So 
far we have shown results for T = Q.le/kB, which is near 
the triple point for spherical Lennard- Jones molecules, 
and roughly 30% above the glass transition temperature 
for chain molecules. This temperature was chosen be- 
cause many airborn hydrocarbons are above their melt- 
ing point at room temperature. Figure ^ shows the vari- 
ation of a and tq with temperature at w = O.OOScr/iij 
for the system of Figs, [l] and |2| The value of a changes 
only weakly with temperature, falling about 20% from 
ksT/e = 0.3 to 1.2. The high end of this temperature 
range is more than twice the glass transition tempera- 
ture for chain molecules and close to the critical tem- 
perature for spherical molecules. 

The weak variation of a with T is consistent with our 
geometrical picture, since the hard-sphere diameter is 
insensitive to temperature. In contrast, the value of tq 
is roughly proportional to temperature over the studied 
range. Equation implies that the friction vanishes 
for P < —To/ a. Of course this can not happen, but 
there is a qualitative change in the type of friction as the 
pressure drops below —To/a. At high pressures we ob- 
serve static friction and a slightly smaller kinetic friction 
that rises logarithmically with velocity. For P < —To /a 
the static friction vanishes, and the kinetic friction rises 
from zero as a power of the velocity. 

This change in frictional behavior reflects a change 
in the state of the film that has been studied previously 
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Figure 5. Variation of (a) a and (b) tq with temperature at 
V = 0.005(T/t£j for the system of Fig. |^. 



in multi-layer systems |40|,^,^,^3|4j,^ . At high pres- 
sures the molecules lock into a glassy state, while at low 
pressures they diffuse freely and act like a viscous lubri- 
cant. Fig. H implies that the pressure required to lock 
molecules into a glassy state rises linearly with T. This 
result is specific to our choice of purely repulsive po- 
tentials (r'^ = Adding an attractive tail to the 
potential shifts the entire temperature dependent curve 
upwards, and gives glassy behavior at lower (even nega- 
tive) pressures. Note that earlier results on multi-layer 
films indicated that r was independent of velocity in 
the glassy state [Q. However, these results covered a 
relatively small range of v and are not inconsistent with 
the logarithmic variation seen in Fig. |^. 

3.3. Relation to static friction 

The changes in kinetic friction with temperature, in- 
teraction potential, geometry and other factors are very 
similar to what we have observed for static friction 
[ p^l|25[] . Indeed the ratio between the value of ak{v) 
at a given velocity and the value of as for static friction 
is relatively constant. For example, over the range of 
temperatures and potential parameters shown in Fig- 
ures § and H, ak{v = 0.005a /tLj) is between 75 and 
85% of a^. This variation is comparable to our statis- 
tical uncertainties. 

Figure^ illustrates the correlation between static and 
kinetic friction as 6 is varied. Only the range 6* = to 
30° is shown, since the values at other 9 are related by 
symmetry. The simulation data in previous figures was 
for 6 — 30° and the commensurate case corresponds to 
6* = 0°. We confirmed that the friction is linear in P 
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Figure 6. Variation of (a) static and (b) kinetic friction with at 
P = 24ea~^. Error bars in both panels are about 5%. For these 
simulations the area of the unit cell is 2880fT-^ , and the number of 
monomers is equal to the number of atoms on the surface of the 
bottom wall. 

for several 9, and show results for a single, relatively 
high pressure (P = 24eCT~^) to minimize the number 
of calculations, and emphasize the contribution from a. 
Data for a monolayer film is shown because it gives the 
largest variation in friction with 9 and thus the best test 
of correlations between Fk and Fs [Q. As 9 changes 
from 2° to 30°, Fk varies by almost a factor of two, 
yet the ratio Fk / Fg stays constant within our statistical 
uncertainties. For v — 0.005(7 /t^j wc find Fk/Fg = 
0.88±0.10, while Fk/Fg = 1.05±0.10 for v = 0.05a/tLj. 

The correlation between Fk and Fg is lost when the 
walls become commensurate. As 9 decreases from 2° to 
0° there is a sharp increase in Fg by about a factor of 
seven, and a drop in Fk by about a factor of two. In 
addition, there is a qualitative change in the velocity 
dependence of the kinetic friction: Fk decreases with 
velocity in the commensurate case, but increases with 
V for incommensurate walls. 

The observed correlation between static and kinetic 
friction is somewhat surprising given that they reflect 
different physical effects. As noted in the introduction, 
the static friction is related to the force needed to pull 
the system out of a local potential energy minimum, 
while the kinetic friction is related to the energy dis- 
sipated during sliding. Nevertheless, the above com- 
parisons of Fk and Fg indicate that the two arise from 
similar processes when the walls are incommensurate. 

3.4- Comparison to the Tomlinson model 

Some insight into the relation between static and ki- 
netic friction is provided by a simple one-dimensional 
ball and spring model called the Tomlinson model 
Jl(J,[lli. This model considers two bare surfaces in 
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contact. The bottom surface is modeled by a sinusoidal 
potential with amplitude Vq and period dnn- Atoms in 
the top surface are treated as independent, overdamped 
oscillators that are coupled to the center of mass Xcm 
by a spring of stiffness k. 

The behavior of the Tomlinson model depends on 
the dimensionless parameter A = 47r^Vo/fcd^„ that mea- 
sures the relative stiffness of the potential and spring. 
When A is small, the spring is so stiff that each atom 
has only one metastable state at any Xcm- In this limit 
atoms move smoothly over the periodic potential from 
the bottom wall. This nearly reversible motion leads to 
little dissipation, and Fk goes to zero linearly with the 
center of mass velocity P Jl0| , pT| . When A is large, each 
atom can get locked in many metastable states at any 
center of mass position. Each atom starts in one such 
state. As Xcm increases, that state eventually becomes 
unstable and the atom pops rapidly to the next state. 
During these pops the system is far from equilibrium 
and a substantial amount of energy is dissipated. The 
sequence of pops and the amount of energy dissipated 
during a given wall displacement AXcm is independent 
of wall velocity in the limit Vcm 0. Since the dis- 
sipated energy equals Fk{v)AXcM, the kinetic friction 
approaches a constant value, -Ffc(O), in this limit. One 
can show that Fk{0) approaches the maximum force ex- 
erted by the periodic potential Fmax = iirVo / dnn times 
the number of atoms A'' as A — >• oo P, p^yTl] ] . 

The static friction in the Tomlinson model depends 
on whether the walls are commensurate or incommen- 
surate In the commensurate case all atoms feel 
the same periodic potential at each Xcm- The maxi- 
mum force that must be overcome as they advance by 
a period is Fs = NF^ax- The ratio Fk/Fg vanishes 
for A < 1 and approaches unity as A ^ oo. In the 
incommensurate case Fg = Q for A < 1. For A > 1, 
the static friction is equal to the low velocity kinetic 
friction. This link between Fg and F^ can be seen by 
calculating Fj. (0) from instantaneous forces rather than 
the dissipated energy. At very low velocities almost 
all of the atoms are stuck in a metastable minimum at 
any instant in time. The force resisting motion, Ffc(O), 
comes from the force exerted by the periodic potential 
on these atoms. One can show that summing this force 
or calculating the energy dissipated during atomic pops 
gives the same answer for Fk (l^jlll ■ For incommensu- 
rate walls the distribution of metastable states and thus 
the lateral force is independent of position. Hence the 
static friction equals the force Ffc(O) needed to advance 
the wall at a very low velocity. 

In most cases two bare surfaces will be incommensu- 
rate and the interactions between them will be weaker 
than their internal interactions. This corresponds to the 
case A < 1 and one expects the static friction and low 
velocity kinetic friction to vanish §|,|,|j^J|,|,0jll) . 
In contrast, the layer of adsorbed molecules included 



in our calculations acts like a very compliant solid with 
A > 1. As expected from the Tomlinson model, this pro- 
duces a finite static friction and a comparable kinetic 
friction. The lateral force is predominantly determined 
by the hard sphere interactions which can be expressed 
as the normal load times the local slope of the surface of 
closest approach (Figure Since the surface is nearly 
independent of load and other parameters, the lateral 
force is proportional to load and Amontons' laws are 
obeyed. 

Wc have confirmed the connection between our re- 
sults and the Tomlinson model by analyzing the motion 
of individual monomers at low velocities. The number 
of particles that move by more than a distance c over 
a time interval of 2tLj is calculated as a function of c. 
This time interval is chosen because it is longer than 
the velocity-velocity autocorrelation time and the typi- 
cal period of oscillations of individual monomers about 
their metastable minima. We focus on the case where 
the film is in a glassy state and static friction is observed 
(P > -To/a). 

In equilibrium {v — 0), the distribution of monomer 
displacements is consistent with thermal motion about 
metastable minima plus extremely slow diffusion of 
monomers between metastable minima. When the top 
wall is advanced at a low velocity, almost all of the 
monomers are locked in local metastable minima dur- 
ing any short time interval. As in the Tomlinson model, 
they provide almost all of the force on the top wall and 
thus determine Fk- Monomers that advance by signif- 
icantly more than thermal displacements are far fewer 
in number and contribute a smaller force per atom to 
Fk. For example, with v = O.OOOScr/iij at ksT/e = 
only 0.3% of the monomers move by more than 0.12(7 in 
2tjL. These monomers contribute less than 0.1% to the 
lateral force on the walls, and the percentage decreases 
as the velocity is lowered. 

When the walls are commensurate, the distribution 
of metastable states depends on the relative positions 
of the walls. The number of atoms that pop between 
metastable states is a periodic function of the wall dis- 
placement. No pops occur when the walls are aligned 
directly above each other, and the total energy of the 
system is minimized in this configuration. The num- 
ber of pops rises rapidly as the walls are moved out of 
alignment, and there are pronounced clusters of pops 
at unfavorable wall positions. The static friction is the 
maximum force exerted on the top wall at any position, 
while the kinetic friction represents an average over all 
wall positions. As a result, the two can be very differ- 
ent and the ratio depends on many factors including the 
length of the chains, number of chains, etc.. 

When the walls are incommensurate, the distribu- 
tion of metastable states is independent of wall position. 
Thus the kinetic and static friction are correlated. The 
analogy to the Tomlinson model would suggest that the 
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kinetic and static friction should be equal in this limit. 
As illustrated above, this ratio depends on the velocity 
at which Fk is calculated and it may also depend on how 
Fs is calculated. Quoted results for Fg were obtained 
by determining the minimum force needed to produce 
a displacement of cr in a time IOO^lj. One might thus 
expect the results to correspond to an effective veloc- 
ity of order O.Ola/tLj, and the static friction is indeed 
between the values of Fk foi^w = 0.005 and O.Oba/tLj- 



Rate-state models (Eqs. 'iA and 3^) also imply that 
the static friction depends on the time spent at zero ve- 
locity before a lateral force is applied. We are currently 
exploring the effect of waiting time on static friction and 
its relation to the velocity dependent kinetic friction. 



4. Summary and Conclusions 

We have examined the kinetic friction of surfaces sep- 
arated by a layer of adsorbed molecules as a function 
of velocity, temperature, interaction potentials and sur- 
face alignment. The results are qualitatively consistent 
with many macroscopic measurements, and can be un- 
derstood by analogy to simple model systems. 

One of our main results is that the kinetic friction 
rises logarithmically with velocity. This dependence is 
commonly observed in experimental systems and is of- 
ten associated with thermal activation ||2^,^,^. Our 
analysis of the dynamics of individual monomers is con- 
sistent with this idea. We found that the kinetic friction 
was given by the force exerted by monomers that were 
trapped in metastable states. A simple argument shows 
that this force will decrease as v decreases. 

The force from the top wall on a monomer in a given 
local potential well is zero at the bottom of the well. 
As it is displaced away from the minimum the force 
rises. When the force reaches its maximum value the 
monomer becomes unstable and pops to a new well. 
Since atoms that exert the largest force are closest to be- 
coming unstable, they are the ones that are most likely 
to be thermally activated. When they hop to a new 
well, they contribute a smaller lateral force. As v de- 
creases the time for thermal activation increases as a /v. 
Thus more and more of the large contributions to Fk are 
lost due to thermal activation. As a result, Fk drops as 
V decreases. 

Another important result is that adsorbed layers nat 
urally lead to a kinetic friction that satisfies Eq. 1 
He et al. had previously established this relation for 
the static friction Thus adsorbed layers provide a 
molecular basis for Bowden and Tabor's explanation of 
Amontons' laws in both contexts. 

Changes in the kinetic friction with interaction po- 
tential and temperature closely parallel changes in the 
static friction. Both Fk and Fg are insensitive to 
parameters that are not controlled in typical experi- 
ments. Examples include wall orientation (9), the num- 



ber and precise size of adsorbed molecules, and even 
the strength of adsorption (cto/). This insensitivity is 
consistent with the fact that variations in the measured 
friction between nominally dry surfaces are relatively 
small 25%). In contrast, simulations and analytic 
studies of bare surfaces show dramatic variations with 
surface orientation . 

The linear relation between shear stress and pressure 
in Eq. |l.l| can be understood from a simple hard sphere 
picture (Fig. ^). The wall atoms create a surface of 
closest approach that adsorbed molecules must climb up 
in order to move laterally. The lateral force and normal 
force are related by the local slope of this surface, and 
a represents some average of the slope. Since e^f has 
little affect on the hard sphere diameter, changing it 
has little affect on a. In contrast, changing the ratio 
f^wf/dnn has a direct effect on both the slope and a 
(Figs. I and I). 

The simple hard sphere model also explains the in- 
sensitivity of a to the number of adsorbed molecules 
per unit surface area p [|l2|,|5|. Increasing the density 
of molecules spreads the normal pressure P and lateral 
stress T over more molecules. However, the ratio be- 
tween them, a, is still determined by the slope of the 
surface of closest approach. This surface is determined 
solely by the solids when the film is less than a mono- 
layer thick. One only expects variations in a when the 
film becomes thicker than a monolayer and can shear 
internally. Preliminary studies of the thickness depen- 
dence do show such effects as the thickness increases to 
two or three layers. 

The close connection between our results for kinetic 
and static friction is typical of experiments on nominally 
dry surfaces and surfaces that are purposely coated with 
thin hydrocarbon films . We discussed how this con- 
nection can be understood from results for the simple 
Tomlinson model. Although the ratio Fk{v)/Fs is rel- 
atively constant for a given velocity, there are notice- 
able changes with v due to the logarithmic variation in 
Fk{v). We expect that the value of Fg will also vary 
logarithmically with the measurement time, and this 
remains an avenue for future research. 

Another important focus of future research will be 
to improve the quantitative comparison between the- 
ory and experiment. The simulations described here 
use simple interaction potentials and surface geome- 
tries that capture qualitative features of the friction due 
to adsorbed layers. More detailed models of molecular 
structure and bonding, as well as surface roughness and 
elasticity will be needed to allow direct comparison to 
experiment. 
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